A Simulation Method to Resolve Hydrodynamic Interactions in Colloidal Dispersions

A new computational method is presented to resolve hydrodynamic interactions acting on solid particles immersed in incompressible host fluids. In this method, boundaries between solid particles and host fluids are replaced with a continuous interface by assuming a smoothed profile. This enabled us to calculate hydrodynamic interactions both efficiently and accurately, without neglecting many-body interactions. The validity of the method was tested by calculating the drag force acting on a single cylindrical rod moving in an incompressible Newtonian fluid. This method was then applied in order to simulate sedimentation process of colloidal dispersions.Comment: 7pages, 7 figure

Numerical methods for multibody systems

This article gives a brief summary of some results obtained by Nasser on modeling and simulation of inequality problems in multibody dynamics. In particular, the augmented Lagrangian method discussed here is applied to a constrained motion problem with impulsive inequality constraints. A fundamental characteristic of the multibody dynamics problem is the lack of global convexity of its Lagrangian. The problem is transformed into a convex analysis problem by localization (piecewise linearization), where the augmented Lagrangian has been successfully used. A model test problem is considered and a set of numerical experiments is presented

A new approach to hyperbolic inverse problems

We present a modification of the BC-method in the inverse hyperbolic problems. The main novelty is the study of the restrictions of the solutions to the characteristic surfaces instead of the fixed time hyperplanes. The main result is that the time-dependent Dirichlet-to-Neumann operator prescribed on a part of the boundary uniquely determines the coefficients of the self-adjoint hyperbolic operator up to a diffeomorphism and a gauge transformation. In this paper we prove the crucial local step. The global step of the proof will be presented in the forthcoming paper.Comment: We corrected the proof of the main Lemma 2.1 by assuming that potentials A(x),V(x) are real value

Exact Controllability of the Time Discrete Wave Equation: A Multiplier Approach

In this paper we summarize our recent results on the exact boundary controllability of a trapezoidal time discrete wave equation in a bounded domain. It is shown that the projection of the solution in an appropriate space in which the high frequencies have been filtered is exactly controllable with uniformly bounded controls (with respect to the time-step). By classical duality arguments, the problem is reduced to a boundary observability inequality for a time-discrete wave equation. Using multiplier techniques the uniform observability property is proved in a class of filtered initial data. The optimality of the filtering parameter is also analyzed

Ségrégation suivant la forme et la taille des grains pour des écoulements dans un tambour tournant

Afin de prendre en compte l'angularité des particules dans les écoulements granulaires, l'IFP a développé un outil numérique : GRAINS3D. Ce dernier repose sur la méthode des éléments discrets et un algorithme de détection de contact développé pour la réalité virtuelle (algorithme Gilbert-Johnson-Keerthi). A l'aide de ce code, nous proposons, dans cette communication, d'étudier et de quantifier l'influence de l'angularité sur les écoulements granulaires. Nous nous intéressons au phénomène de ségrégation pour des massifs granulaires morphologiquement homogènes qui contiennent des particules de taille différente et des mélanges binaires de formes. Les résultats obtenus sont confrontés à ceux de la littérature. Dans le cas de formes plus complexes, une analyse du comportement des particules en fonction de leur convexité est proposée

DNS of vertical plane channel flow with finite-size particles: Voronoi analysis, acceleration statistics and particle-conditioned averaging

We have performed a direct numerical simulation of dilute turbulent particulate flow in a vertical plane channel, fully resolving the phase interfaces. The flow conditions are the same as those in the main case of "Uhlmann, M., Phys. Fluids, vol. 20, 2008, 053305", with the exception of the computational domain length which has been doubled in the present study. The statistics of flow and particle motion are not significantly altered by the elongation of the domain. The large-scale columnar-like structures which had previously been identified do persist and they are still only marginally decorrelated in the prolonged domain. Voronoi analysis of the spatial particle distribution shows that the state of the dispersed phase can be characterized as slightly more ordered than random tending towards a homogeneous spatial distribution. It is also found that the p.d.f.'s of Lagrangian particle accelerations for wall-normal and spanwise directions follow a lognormal distribution as observed in previous experiments of homogeneous flows. The streamwise component deviates from this law presenting significant skewness. Finally, a statistical analysis of the flow in the near field around the particles reveals that particle wakes present two regions, a near wake where the velocity deficit decays as 1/x and a far wake with a decay of approximately 1/(x*x).Comment: accepted for publication in Int. J. Multiphase Flo

On discretization in time in simulations of particulate flows

We propose a time discretization scheme for a class of ordinary differential equations arising in simulations of fluid/particle flows. The scheme is intended to work robustly in the lubrication regime when the distance between two particles immersed in the fluid or between a particle and the wall tends to zero. The idea consists in introducing a small threshold for the particle-wall distance below which the real trajectory of the particle is replaced by an approximated one where the distance is kept equal to the threshold value. The error of this approximation is estimated both theoretically and by numerical experiments. Our time marching scheme can be easily incorporated into a full simulation method where the velocity of the fluid is obtained by a numerical solution to Stokes or Navier-Stokes equations. We also provide a derivation of the asymptotic expansion for the lubrication force (used in our numerical experiments) acting on a disk immersed in a Newtonian fluid and approaching the wall. The method of this derivation is new and can be easily adapted to other cases

Numerical analysis and simulation of the dynamics of mountain glaciers

In this chapter, we analyze and approximate a nonlinear stationary Stokes problem that describes the motion of glacier ice. The existence and uniqueness of solutions are proved and an a priori error estimate for the finite element approximation is found. In a second time, we combine the Stokes problem with a transport equation for the volume fraction of ice, which describes the time evolution of a glacier. The accumulation due to snow precipitation and melting are accounted for in the source term of the transport equation. A decoupling algorithm allows the diffusion and the advection problems to be solved using a two-grids method. As an illustration, we simulate the evolution of Aletsch glacier, Switzerland, over the 21st century by using realistic climatic conditions

Continuum viscoplastic simulation of a granular column collapse on large slopes: μ(I) rheology and lateral wall effects

We simulate here dry granular flows resulting from the collapse of granular columns on an inclined channel (up to 22°) and compare precisely the results with laboratory experiments. Incompressibility is assumed despite the dilatancy observed in the experiments (up to 10%). The 2-D model is based on the so-called μ(I) rheology that induces a Drucker-Prager yield stress and a variable viscosity. A nonlinear Coulomb friction term, representing the friction on the lateral walls of the channel, is added to the model. We demonstrate that this term is crucial to accurately reproduce granular collapses on slopes ≳10°, whereas it remains of little effect on the horizontal slope. Quantitative comparison between the experimental and numerical changes with time of the thickness profiles and front velocity makes it possible to strongly constrain the rheology. In particular, we show that the use of a variable or a constant viscosity does not change significantly the results provided that these viscosities are of the same order. However, only a fine tuning of the constant viscosity (η=1 Pa s) makes it possible to predict the slow propagation phase observed experimentally at large slopes. Finally, we observed that small-scale instabilities develop when refining the mesh (also called ill-posed behavior, characterized in the work of Barker et al. [“Well-posed and ill-posed behaviour of the μ(I)-rheology for granular flow,” J. Fluid Mech. 779, 794–818 (2015)] and in the present work) associated with the mechanical model. The velocity field becomes stratified and the bands of high velocity gradient appear. These model instabilities are not avoided by using variable viscosity models such as the μ(I) rheology. However we show that the velocity range, the static-flowing transition, and the thickness profiles are almost not affected by them